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olution of the Radiative Transfer in the 
reflective domain for non absorbing 
atmosphere and lambertian ground 
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(= albedo for lambertian) 
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Perfect Lambertian Reflector 

/ Radiance of the Perfect Lambertian Reflector 

n 2 n s' 

[ f RPLF{0 S , 0, f) cos(#) sin (0)d0df = E s cos(0 s ) 
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Isotropic 

radiation 


P Perfect Lambertian reflector (&s > &v >v) = 1 
P Lambertian reflector (P S ,0 V ,</>) = P 
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Simple Radiative Transfer 

Equation 


ground 


Atmospheric Correction of Earth Observation Data for Environmental Monitoring: Theory and Best Practises 




SRTE (cont.) 
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SRTE (cont.) 
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SRTE (cont.) 
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SRTE (cont.) 
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SRTE 1 interaction (cont.) 


ground 



E 0 _ _ TJOJE^ _ 

E, E, E, 


T (6 V ) P ground 1 (6 s) 


Atmosphefi 



rrect /orfW Earth for EVfv 


vironmental Monitoring: Theory and Best Practises 




SRTE 2 interactions 
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SRTE Multiple Interactions 
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STRE for non absorbing atmosphere 

and lambertian ground 


Atmospheric reflectance Ground reflectance 

(= albedo for lambertian) 
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The composition of the 

atmosphere 


Permanent Constituents 

Variable constituents 

Constituent 

% by 

Constituent 

% by volume 


volume 



Nitrogen (N 2 ) 

78.084 

Water Vapor (H 2 0) 

0.04 

Oxygen (0 2 ) 

20.948 

Ozone (0 3 ) 

12 x 10" 

Argon (Ar) 

0.934 

Sulfur dioxide (S0 2 ) b 

0.001 x 10" 

Carbon dioxide (C0 2 ) 

0.033 

Nitrogen dioxide (N0 2 ) 

0.001 x 10" 

Neon (Ne) 

18. 18 x 10" 

Ammonia (NH 3 ) 

0.001 x 10" 

Helium (He) 

5.24 x 10" 

Nitric oxide (NO) 

0.0005 x 10" 

Krypton (Kr) 

1. 14 x 10" 

Hydrogen sulfide (H 2 S) 

0.00005 x 10" 

Xenon (Xe) 

0.089 x 10" 

Nitric acid vapor 

trace 

Hydrogen (H 2 ) 

0.5 x 10" 



Methane (CH 4 ) 

1.5x10" 



Nitrous Oxide (N 2 0) 

0.27 x 10" 



Carbon Monoxide (CO) 

0.19x 10" 
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Gaseous Absorption (H 2 0) 
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Modified SRTE to account for 

absorption 

In case of a pure molecular atmosphere (no aerosol) we can write: 


Pa Pp (0,A,<f>)=Tg 


othergases k 


r 


Pal J0 s , 0 r ,t) + Tg H i°{m,U^ 0 )T at (e s )T at (d v ) 


gaz> 


P ground 


h 2 o . 


^ ^ a tmP ground 


m is the air mass = l/cos(0 s )+l/ cos(0 v ) 
U is the gaz concentration 
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Final SRTE approximation 
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Spectral response 


Water vapor effect for different 
sensors in the near infrared 


Spectral response of the near infrared band 
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Scattering angle , 

• The scattering angle, ©, is the relative 
angle between the incident and the 
scattered radiation 
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Phase function 


The phase function, P(©), describe the 
distribution of scattered radiation for one 
or an set of particles. It is normalized 
such as: 

\\P(@)dcQ = 4n 


0 0 


2nn 


n 


since 


J| P(0)sin(6)d6d{ it =2;rJ P(0)sin(0)d0 


0 0 
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we have J P(fy sin(0)dQ = 2 
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Rayleigh/molecular scattering 

1/4 

• Rayleigh or molecular scattering refers 
to scattering by atmospheric gases, in 


that case: 


P(0) = — (l + cos 2 (© )) 
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Rayleigh/molecular scattering 2/4 

• The concentration in scatterer is better described by 
the efficiency they scatter at a certain wavelength or 
the proportion of direct transmission which is related 
to the spectral optical thickness x(k) 

E„(X) 

▼ 

E,(X)/E 0 (X)=e- T ( > -> 


.. E t (X) 

• For Rayleigh x(X) is proportional to At 4 and for standart 
pressure is ~ 0.235 at 0.45 pm 
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Rayleigh/molecular scattering 

3/4 

The rayleigh reflectance, p R , could be 
crudely approximated by: 



4p(®) 

4cos(6> )cos(<9 v ) 
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Rayleigh/molecular scattering 4/4 


Compute the reflectance of the sky 
(assumed clear no aerosol) at solar 
noon at 45degree latitude at vernal 
equinox looking straight up at 0.45pm, 

0.55pm, 0.65pm 

p° R A5 ' on (4 5 o,O o ,0) = 


0.235/ > (45°) 


4 cos(45 °) cos( 0°) 


0.235x0.75(1 + cos 2 (45)) 0.265 

4 cos(45°) cos(0°) 2.8 

p° R 55/Jm = 0.1^(0.45/0.55) 4 ~ 0. 1x0.4= 0.04 

p° R 65fJm = 0.04x(0.55/0.65) 4 ~ 0.04x0.5 =0.02 
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Aerosol scattering 1/5 

• aerosol scattering refers to scattering by particles m suspension in the 
atmosphere (not molecules). The MIE scattering theory could be 
applied to compute the aerosol phase function and spectral optical 
depth, based on size distribution, real and imaginary index. 
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Aerosol scattering 2/5 


Continental aerosol phase function 

itf ================= 



20 40 60 80 100 120 140 160 180 

Scattering angleS 

0.400 0.488 0.515 0.550 0.694 0.860 2.250 3.750 
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Aerosol scattering 3/5 



single scattering albedo 

(0.2- 1.0) to account for absorbing 

particles 



4cos(0 4 )cos(< 9 v ) 
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Aerosol scattering 4/5 

Continental aerosol optical thickness spectral variation 
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Aerosol scattering 5/5 

Continental aerosol single scattering albedo spectral variation 
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Atmospheric effect: Vegetation 1/3 


Surface reflectance (Kuusk’s model) 
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Atmospheric effect: Vegetation 2/3 


Top of the atmosphere reflectance (Kuusk/s model) 


1.100 

1.000 

% 0.900 

S 0.800 

? 

0.700 

0.501600 

0.447 

0 . 3^ 00 

0.337 

0.283 400 

0.228 



No absorption, Continental aerosol 

Atmospheric Correction of Earth Observation Data for Environmental Monitoring: Theory and Best Practises 



Atmospheric effect: Vegetation 3/3 


Top of the atmosphere reflectance (Kuusk/s model) 
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Atmospheric effect: Ocean 1/2 

Surface reflectance (Ocean model) 
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(^ 0 , 


Atmospheric effect: Ocean 2/2 


Top of the atmosphere reflectance (Ocean model) 
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Perfect Lambertian Reflector 


n 2n 


f f RPLF(0 s , 0, </>) cos(O) sin( 6)d0d(j) = E s cos(O s ) 

0 0 


9 S 
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Isotropic 
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Different Types 



Specular reflector (mirror) 



Reflectors 



diffuse reflector (lambertian) 



Nearly Specular reflector (water) 




nearly diffuse reflector 
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Gray level temperature image 


MODIS 
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MODIS data illustrating the hot-spot over dense vegetation 
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Scattering Angle (°) 



BRDF atmosphere coupling correction 


Lambertian infinite target approximation 



BDRF atmosphere coupling approximation 




app r atm V s *> v 

r Tl Vvf (a \ ~^ i 
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Adjacency effect correction 


Lambertian infinite target approximation 
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T atm(A) T atm(A) 



Pi 

^ atm P i 



Pi = 


PP T '_^_ + Ctm(A)Pe 
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Adjacency effect correction (testing) 



Synthetic data set for surface reflectance in the blue (a), green (b) and red (c), and in 
RGB (d) corresponding to bright soil (yellow squares) and dense vegetation (green 
squares). 



(a) 


(b) (c) (d) 


Typical atmospheric effect on the synthetic surface reflectance shown above 
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Adjacency effect correction (testing) 



Reflectance’s observed over a horizontal transect on the checkerboard. The red bars are the “true” surface 
reflectance, the blue bars correspond to the top of the atmosphere signal including adjacency effect. The green 
bars correspond to the corrected data using the infinite target assumption. The open square correspond to the 
data corrected for the adjacency effect using the operational method developed. 
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Adjacency effect correction (validation) 
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